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ATTRACTION AND STABILITY OF NONLINEAR ODE'S USING 
CONTINUOUS PIECEWISE LINEAR APPROXIMATIONS 
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Abstract. In this paper, several results concerning attraction and asymptotic 
stability in the large of nonlinear ordinary differential equations are presented. 
The main result is very simple to apply yielding a sufficient condition under 
which the equilibrium point (assuming a unique equilibrium) is attractive and 
also provides a variety of options among them the classical linearization and 
other existing results are special cases of the this main theorem in this paper 
including and extension of the well known Markus-Yamabe conjecture. 

Several application examples are presented in order to analyze the advan- 
tages and drawbacks of the proposed result and to compare such results with 
successful existing techniques for analysis available in the literature nowadays. 



1. Introduction 

Stability of nonlinear continuous system of Ordinary Differential Equations (ODE 
system) is a very important subject either in control design or pure theoretical anal- 
ysis. It is always desirable to have systematic, necessary and sufficient conditions 
able to predict any class of stability where Lyapunov stability is the most common 
definition. 

Among the available techniques nowadays, we can count mainly on Lyapunov 
functions as the most general method whereas many particular cases including all 
the available analysis for linear systems and local analysis also received considerable 
attention, see [T] for a general reference in Lyapunov classical theory and [37] j US] 
for successful applications. The well known Lyapunov theorem suffers the lack of 
a systematic counterpart making it impossible to determine if such a Lyapunov 
function even exists for general cases. For instance, the example by Krasowski in 
[3], pp. 32, |24| and [5] shows that sometimes such a continuous Lyapunov functions 
does not exists. 

Several tools to construct Lyapunov functions were derived in the literature, for 
instance, the special case in [33] provides a methodology for polynomial systems 
with constraints, on the other hand, the paper in |18| was the first constructive 
converse Lyapunov method, in fact it was proved there that if an equilibrium point 
posses exponential stability behavior, then a certain feasible linear programming 
problem yields a Lyapunov function for such a system. Later, the work in jl9| 
improved the result proving a similar theorem for asymptotic stability instead of 
exponential one. Moreover, the important reference [35] shows that the analysis of 
Lyapunov functions using polynomial approximations is topic of interest nowadays. 
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Alternative techniques to determine if a given system is stable (asymptotic or 
attractive) are the method conjectured by Aizerman (see |38) and [5], pp. 86, 
Theorem 2) and the result in p5(J; where is required to be concave. Also the 

Corollary 2.43 in the book [8' is an alternative to the Lyapunov's method, in fact 
this result only requires to check the stability of the jacobian of a given nonlinear 
autonomous system at the origin (equilibrium point), the theorem also proves that 
if there exists a negative upper bound for all the real part of the eigenvalues of 
the Jacobian matrix, then there exists a region surrounding such an upper bound 
where the system presents exponential behavior 

Particular cases of ODE's can be successfully treated specially for low dimen- 
sions; namely and 3?^ or even discontinuous systems with linear components 
which can be classified under the names of Hybrid Linear, Piecewise Linear, Switch- 
ing Linear, etc. For these cases, several techniques have been developed; such as. 
Common Quadratic Lyapunov Function, Multiple Lyapunov Functions, Poinare 
Maps, Lie Algebra among many others (see for instance [21], [22], [10], [27], [26] y 
[28]). 

Besides the interest in stability or attraction, and as it is well known, it is neither 
simple nor systematic to determine a region where any initial condition drives the 
system to an stable equilibrium point (domain of attraction). The most general 
analytic framework to find the domain of attraction requires to solve a Halmiton- 
Jacobi type partial differential equation (see [3 ), in this regard the paper in [36] 
provides such a region by defining Lyapunov (like functions which allows to prove 
that the trajectories never leaves the region), whereas the classical Lyapunov func- 
tions only prove the existence of such a region without providing it. Finally notice 
that according the section Related Work in that paper, only two methodologies are 
systematic nowadays: [33] and [34] . 

In particular the method presented in f6^ improves the early algorithm in [361 re- 
moving some computational drawbacks, however, the algorithm starts with a given 
Lyapunov function, otherwise, if such a function is not provided then a linearization 
at the origin along with a quadratic form is used to run the method, it turns out 
that for certain cases (called critical cases) the linearization exhibits zero or pure 
complex conjugate eigenvalues making the construction of the quadratic Lyapunov 
function more difficult with the necessity of some extra machinery [1 3J . 

The present paper presents two type of sufficient conditions, namely, a sufficient 
condition for asymptotic stability in nonlinear autonomous ODE's which makes 
uses of a constant negative definite matrix to produce two outcomes: 

• A quadratic Lyapuonv function for the nonlinear ODE being analyzed. 

• A subset of the domain of attraction. 

the other sufficient condition is a corollary to the previous main result, in fact 
instead of using a constant stable matrix, it uses the jacobian of the nonlinear 
vector field as a function of the state-space variables, in this case the only outcome 
is a subset of the domain of attraction. 

The results in this paper are built on the basis of the error bounds derived in 
[16] and |14| . where given a Nonlinear autonomous ODE, then an approximation 
of that vector field using a Continuous Piecewise Linear (CPWL) basis is consid- 
ered providing dynamic bounds for an error defined as the difference between the 
trajectories of the Nonlinear ODE and its CPWL approximation. 



It turns out that the drawback is to determine such a region. 
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This dynamic bounds (error bounds) posses the same dynamic as the CPWL 
approximation ODE, in this way, stabihty issues for the CPWL ODE are also 
shared by the error bounds. ReaUze that once the stabihty for the CPWL ODE is 
proved, then stability for the dynamic bounds is also ensured then the trajectories 
of the Nonlinear ODE tends to a constant or in other words also shares stability 
issues. 

It is worth to mention the works [12] and [11] . where interval analysis it is used 
in order to derive a method, both for asymptotic stability and the determination 
of some subset of the attraction domain. In fact, in |11| . an algorithm is provided 
to check asymptotic stability for a nonlinear ODE using successive linearization 
around the equilibrium at the origin. However, the stability proof is based on a 
quadratic Lyapunov function using the linearization aforementioned, which is not 
ensuring the algorithm ends in finite time for general nonlinear cases. 

The work [S] is very interesting and a new result, in fact, the authors of that paper 
derive necessary and sufficient conditions under which a given nonlinear ODE is 
asymptotically stable. In this way, the results in the present paper can be compared 
to the results in 9 to notice that the present paper posses several results more 
general and also including the ones in that paper, however the tools in 9^ can 
be utilized together with the results in the present paper, then allowing a more 
extensive analysis for nonlinear systems. 

This paper is organized as follows: Section [2] introduces formally the problem 
considered; Section lOl proves a sufficient stability condition for CPWL ODE's while 
Section [3] provides a sufficient condition for stability or attraction in a nonlinear 
ODE using CPWL theory. Finally, Section |4] shows some examples of application. 

2. The Problem Considered 

The goal of this paper is to derive sufficient conditions for attraction and as- 
ymptotic stability as well as quadratic Lyapunov functions for nonlinear ODE's: 
x{t) = f{x). The tool for such an achievement is an extension of the error bounds 
in [12] with the formalism in [T3] and [12] . 

Thus, according to [I^, if a given ODE x{t) — f{x) produces a CPWL approxi- 
mation, such that this approximation posses an asymptotically stable equilibrium 
point, then the error bounds in the mentioned reference, also posses an asymptot- 
ically stable equilibrium point, showing that the properties of stability are shared 
by the given ODE and its approximation. 

A correction term in the development in [16j is added in this paper to take into 
account the possibility for the trajectories of x{t) = f{x) and its CPWL approxi- 
mation to run different simplices during some interval of time. 

In this way, the scenario is as follows: 

x{t) = f{x) Given ODE 

{t) = AC^) -xcpwl^ SC^) CPWL Approx. 
X, xcpwL e W 

Recalling that the goal of the paper is to derive stability conditions of nonlinear 
ODE's using its CPWL approximation in the limiting caso of the grid size tending 
to zercO, this is only possible if the trajectories of both: nonlinear and CPWL are 



^See | 23| for a definition of grid size. 
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close each other. To make precise what close each other means, the concept of error 
bounds used in [16^ and [15, will be utilized. 

In this way, the formalism in [TS] yields the following error bounds: 



'x{t) = f{x) Given ODE 

icPWLit) = AC^) ■XCPWL+ SC^) CPWL Approx. 
I mzn{i,2} {E*'it),E*^t)} < E^'^\t) < max^,^^} {E*\t) , E*^t)} 
^ ' ' I = A^'^ ■ E*^ + f '=(<) - A 

E*-^{t) •S*2 + f'=(<) + A 

i = 1, ..,71 

where e 3fJ"^i, A e sft"xi,^(*'=) = (^(fe) - A(^)).xcpwl + (s('=) -BW) 

and E^^\t) — x{t) — xcPWL{t) is the error vector. 

Notice that extra term S}-^^^ (t) takes into account for the possibility of the tra- 
jectories x{t) running the simplex j*'' and xcpwl running the simplex fc*'*. 

On the other hand, to reach the goal in this paper, it is necessary to make 
sure that bounds in mm{i^2} {E*^ [t] , E*"^ {t)] and ■max{i^2} {E*^ {t) , E*'^ {t)] 

do not grows to infinity. Clearly, if the dynamics of such a bounds: E*-^''^{t) = 
j[M . ^*i.2 _|_ ^ik _|_ ^ j^j.g attractive to one equilibrium point, then in each simplex 
the stability issues of the CPWL ODE is shared by the nonlinear one. 

Since the error bounds in ()2.ip can be rewritten in matrix block form: 



(2.2) 



'E*^^^{t)' 




■^(0 


- A{i) 

















AC) 



±A 




The theorem in next subsection asserts that a CPWL ODE is attractive to one 
of its equilibrium points once all the matrices A*-*^-* are definite negative. In this 
context, (|2.2p will conduct to the desired goal in this paper in the view of the well 
known Theorem for block matrices: the eigenvalues of matrix A^''^ is the junction 
of the eigenvalues of matrices {A'-'^-', A^*-'}. 

2.1. A Sufficient Condition for Asymptotic stability in CPWL ODE's. 

In this section a theorem concerning attraction for CPWL ODE's is recalled. In 
fact, the theorem was proved in [14 and it is important because it provides a tool 
which will be invariant when the grid size of the CPWL approximation simplicial 
partition tend to zero (as it will be shown in section [3]). 



Theorem 1. Given a CPWL ODE: 



(0 



defined in an set 



57 which posses no equilibrium points over the frontier^. If, moreover, for each 
simplex i*^, A^"^^ is negative definite and fl is an invariant se0, then the CPWL 
ODE is attractive to at least one of its equilibrium points wiht the initial conditions: 

x{0) e n. 



This means that the equiUbriums of each individual hnear system in each simplex are outside 
the borders. 

^See | 32| . pp. 73 and pp. 76 for a further reading. 
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Whit this Theorem and considering the grid size of a CPWL ODE approximating 
a given nonhnear ODE, tending to zero (the CPWL approximation tend to the 
nonhnear vector field, '23') , it is going to be proved a theorem providing sufficient 
conditions for attraction and stability in a systematic way. 

3. The Main Result: A Sufficient Condition for Stability of 

Nonlinear ODE's 

Up to now, conditions for a CPWL ODE being asymptotically stable were de- 
rived, in particular when these ODE's come from a CPWL approximation of a 
nonlinear vector field, the quality of the approximation improves as the grid size 
tend to zero (see |23]). 

In order to apply this results to predict stability /attraction in Nonlinear ODE's, 
two main research streams can be followecfl: 

• Given an ODE: x{t) = f{x) and a CPWL approximation of it, analyzing 
the stability of each matrix A^*^ (this is the research line presented in |22j). 

• Given an ODE: x{t) — f{x), utilizing the CPWL theoretical tools (like 
the ones in [23 ), to derive conditions under which there exists a CPWL 
approximation (without explicit calculation of it), such that each matrix 
A^*' is Hurwitz, then ensuring stability or attraction in the view of Theorem 

HI 

Following the second research line above, it is clear that the sufficient conditions 
in Theorem [U require the following: 

• To derive conditions using f{x) such that there exists a CPWL approx- 
imation for certain degree of accuracy where each matrix A^^'> is definite 
negative. 

• To find an invariant set for x{t) = f{x) 

In order to provide a smooth introduction for the main result in this paper it is 
needed some extra machinery. In this way, the first condition above, requires the 
following Lemma: 

Lemma 2 (Stability Conditions for Linear Systems). Given a Linear ODE x{t) — 
A ■ X + B such that: 

-F,k+A,k + B, \<X* 
Bj \< A* 
j,fc = 1,2, . . . ,n 
F' ■ P + P ■ F = -I 
R^\\P\ -(A* + A*) + (A* + A*)'- I P I 
Xji < 1 The biggest eigenvalue of R 
I • I means the absolute value of each element 
Then, A is negative definite-see [4^ for details on negative definiteness. 

Proof. The proof is in the Appendix. □ 



^From now on stability will mean stability in the sense of Lyapunov (see 1321 . pp. 76-82 for 
definitions of Lyapunov stability) . 
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It is important to remember that a positive definite matrix P and symmetric 
such that F' ■ P + P ■ F = —I, always exists if and only if F is Hmwitz (see [25] . 
pp. 22). 

On the other hand, an important corollary which will be useful in what follows is 
the fact that V{x) — x' ■ P ■ x constitutes a Lyapunov function for x{t) = A ■ x + B: 

Corollary 3. Under the hypothesis of Lemma\^ '^{^) — x' ■ P ■ x is a Lyapunov 
function for x(t) — A ■ x + B . 

Proof. The proof is in the Appendix. □ 



Using these results for pure Linear systems and taking into account that all the 
development in this paper is supported by CPWL approximations, which are in 
essence a collection of linear systems, it is possible to prove the main result: 

Theorem 4 (Main Result). Let x{t) = f{x), , /(O) = 0, f e a Nonlinear 
ODE such that f{x) is Lipschitz continuous in a set f2. // given three matrices 

{F,A*,A*}.- 



F^[Fi F2 ... F„] 

A* -[5 5 ... a;] 
x*^[xi A* ... a;] 

F' ■ P + P ■ F = -L 

i?= I P I -(A* + A*) + (A* + A*)'- I P I 

Xr < 1 The biggest eigenvalue of R 

I • I means the absolute value of each element 



f2 = {a 



dfix) 



X \< mini 



X* 



dfix) 
dx 



F |< (A* - A*) 



where Xj is the maximum error between f{x) and the CPWL approximation in 
the spirit of ^ and e^ = [0 ... . . . 0]'. 

i^^position 

Then, x{t) = f{x) is attractive to the origin in Q. Moreover, n is a domain of 
attraction. 



Proof. The proof is in the Appendix. 



□ 



This Theorem provides sufficient conditions for attraction in a region ft, so in 
order to extend this result to asymptotic stability, it will necessary the following 
Lemma: 

Lemma 5 (Lyapunov Functions). Given an ODE: x{t) — f{x), f e f Lips- 
chitz continuous with /(O) ~ 0. If given three matrices {F, X* ,X*}: 
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F ■ 

X* 
A* 



--[XI 
.[XI 



X% 



Fn] 
K] 

Xl\ 



F' -P + P- F ■ 



R = 



then V{x) 
follows: 



P I -(A* + A*) + (A* + A*)'- I P 
Xr < 1 The biggest eigenvalue of R 

X ax 

x' ■ P ■ X is a Lyapunov function for f{x) for all x G Q,with as 



{x : — ■ X \< mini=i, 

ox 



dfjx) 
dx 



F |< (A* - A*) 



Proof. The proof is in the Appendix. 



□ 



This Lemma defines a set f2 where a Lyapunov function V{x) = x' ■ P ■ x is 
available. However, as it is well known, the set il could be a subset of the truly 
domain of attraction. In this way, to find the complete domain of attraction, the 
optimization algorithm in [6] could be used, with the Lyapunov function V{x) = 
x' -P -x. 

Clearly, Theorem H] and Lemma [5] prove trivially the following result: 

Theorem 6 (Sufficient Conditions for Asymptotic Stability). Under the conditions 
on Lemma\^ x[t) = f{x) is asymptotically stable in the set: 



11 ■ X \< mmi=i^,,,^n Aj} 



In this way, the first possible choice for F which always yields a non-empty region 
n containing the equilibrium at the origin, is = ^^^[x=o- This is not more than 
the classic linearization at the origin but this time, with a known validity region. 

A second choice for is: F — providing Sufficient and systematic condi- 

tions under which x(t) = f{x) is attractive to the origin: 

Corollary 7 (Sufficient Conditions for Attraction). Let x{t) = f{x), , /(O) = 
0, f ^ a Nonlinear ODE such that f{x) is Lipschitz continuous in a region fl 
given by: 



= {x : Xb.{x) < 1} 

df(x) 

where R = 2-(| P | • | — ^ • a; | -[1 1 ... 1] +H{xy) and Xr is the biggest eigen- 

^ V ' 

H{x) 

value of R providing P ■ + • P = ^I- Then, x{t) — f(x) is attractive to 

the origin for all x € Q. 

On the other hand is interesting to compare our test of stability in Theorem [3] 
with Theorem 7.1 in [52], pp. 115-116, where the following result is presented.: 
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Theorem 8. Let z(t) he a piecewise trajectory of the system x{t) = f(x) and 
assume that: 

If there exists numbers yi > and symmetric matrices U, W" with non-negative 

entries, and a symmetric matrix T such that Pi = F*' ■T-F'-^^ and Pi = f''^^ -T-F^ 
E TT Fi satisfy: 

{— 2 • • 7^*^ • / > A^*) • F^*) + • A'*^ + £^(*^' • C/(*^ • 
for i in the first simplex and: 



-2.£( 



r(«) 



■E 



■ r -I > A''' ■ p"' + ppw • A 

^E"^' ■ -E* < P(*) < 7W ■ / 
for i in the next simplex, then x{t) tends to zero exponentially. 

Notice that in order to use this theorem and as is pointed out in [22] , we need to 
run a CPWL approximation numerically and then extracting from that the matrices 
to conclude stability for our given nonlinear ODE, the point is how to 
ensure a correct simplicial division to get precise conclusions. In that sense our 
result is much more stronger since we only need the evaluation of the matrix F 
and the Lipschitz coefficients if an estimation of the domain of attraction has to be 
determined, without the necessity of matrices A^-^\ i?^*^. 

We leave this section with a example of application of Corollary [7] and Theorem 
|4] utilizing the jacobian at the origino 

Example 

The well known system of Hopf in normal form is: 

\x = fi{x,y) = a ■ X - y + X ■ {x'^ + y"^) 
\y = f2{x,y) ^ X + a ■ y + y ■ {x"^ + y^) 
Calculating the jacobian matrix: 



dm 

dx 



a + 3 ■ + 
1 + 2- X -y 



-1 + 2 -x-y 

a + x^ +y^ 



then, the Jacobian at the origin leads: 

dfjx) ^ 
dx 



The eigenvalues are a± \/( — 1), which indicates the necessity for a < 0. In this 
way, the conditions on Theorem 2] yields: 



(3.1) 



dfjx) 
dx 



F \< X* - X* ^ 



3-x^ 



y 



2-x-y 



2-x-y 3 • y2 



■ x\< mini^i^2 A* ^ 



<X* -X* 



{a ■ X - y) + 3 ■ X ■ {x^ + y'^) 
{a ■ y + x) + 3 ■ y ■ (x^ + j/^) 



< mini=i^2 A* 



^More examples are given later in Section |4] 
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On the other hand, if A*, A* are such that: 

(X* = [XI XI] 
\X* = [XI XI] 

where A* G and A J € Then, i? is a matrix with rank 1, so the 

condition A^j < 1, yields: 

{A/j = trace{R) < 1 
i?=i-[(A*+A*) + (A*+A*)'] 

where the operator trace{-) is the sum of the diagonal elements. The condition 
for Xr <1 leads: 



(3.2) 



where A* = 



\r<I^Y (Ati + A^a + Ki + A22) < 1 



^11 ^12 
^21 ^^22 



and A* = 



^21 '^22 



Finally, the set of conditions to satisfy arise from equations p.ip and p.2 



(3.3) 



3 -2:2+^2 |<AJi-AJi 
x2+3-y2 |< A^2^-A^2 
2-x-y\<Xl^-Xl^ 

{a ■ X - y) + ?> ■ X ■ {x^ + y'^) |< mini=i^2 

(a • y + x) + 3 • y • (x^ + |< rnin^^i^2 >^2z 

-^(Ati+A^2 + Ati+A52) <1 



Noticing that the condition on Xu is limiting the values for {X*j, X*j} and {A22, Aj^, AI2, A21} 
are not present there, then it is possible to choose any value for them. If in partic- 
ular they are such that: 



-00 

■'21 ^ 

In this way, conditions in equation (13.31) . lead: 




3-.T2+y2 |< A*^ 



3-y2 |<A: 



22 



5i 

A 



22 



-0^08,0. 



< 1 
2- A? 



leading a 



\ia-x-y) + S-x-{x^ +y^) |< A^i 
I (a-y + x)+3-y^(.T2 + y2) |< A^a 

. ~ a ■ (-^11 + ^^22 + '^ll + ^^22) 

One possibility is to choose 2 • A^^ ~ A^j ~ 
domain of attraction a; e [—0.08, 0.08], ye^ 

Next section is providing a set of examples to show the applicability of both: 
Theorem [4] and Corollary [71 together with some discussion about the accuracy 
when the domain of attraction obtained has to be used as an estimation to the 
truly global one. 
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4. Examples 

This section is devoted to show the apphcabihty of the main theorem (Theo- 
rem 2]), as well as the its corollary. They together allow, in some cases, a precise 
estimation of the domain of attraction while asymptotic stability is established. 

Example 1 (An example with finite domain of attraction). This first example is 
taken from the book [3], pp. 53: 



xi{t) 

X2{t) 



-xi + 2 ■ xf ■ X2 

-X2 



For this system is known that the domain of attraction is given by xi ■X2 < 1, then 
this is a boundary for the estimation with Theorem 21 Calculating the Jacobian: 



(4.1) 



dfjx) 
dx 



-1 + 4- xi ■ X2 2 • xf 
-1 



Then, calculating at the origin: 

dm^\-i 

dx [ -1 
This Jacobian is clearly Hurwitz with matrix P as follows: 



1 
1 



On the other hand, the conditions on Theorem |4] leads: 
(4.2) 



dfjx) 
dx 



dfjx) 
dx 



F\< X* -\* ^ 



< A* - A* 



• a; I < mini. 



A ■ xi ■ X2 2 ■ x\ 


Xi + A ■ x\ ■ X2 + 2 ■ x\ ■ X2 
-X2 



< mini 



Analyzing the condition for Xr: 

i?=i-[(A*+A*) + (A* + A*)'] 
On the other hand, if A*, A* are such that: 

(x* = [XI XI] 
\X* = [XI XI] 

where A* £ 5R^^^ and A* G 5R^^^. Then, i? is a matrix with rank 1, so: 

Xji ~ trace{R) < 1 
In this way, the second condition to satisfy arises: 



(4.3) 



Afl < 1 ^ ^ • {Xll + ^22 + Kl + A22) < 1 
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Figure 1. Domain of attraction comparison using Theorem 



where A* 



and KB 



An 
A 



21 



^12 
^22. 



and A* 



A, 



21 



212 
^22 



In summary, from equations 



(4.4) 



4:-xi-X2\< X*ii -X*ii 



2 ■ x\ < A*2 



Al2 

I —Xi + 4: ■ xl ■ X2 + 2 ■ xf ■ X2 \< mirii^ 
I X2 \< mini^i^2 X 



2i 



1 

^. 2 



(An + ^22 + Ki + ^22) < 1 
If in particular, it is chosen: 



I XI2 +00 
[A^i ^ +00 

Then, the conditions in equation (|4.4I) . yields: 



4- I xi • X2 \< All - All 

I a:i I • I -1 + 4 • a;2 + 2 • a;i • 0:2 |< All 

I X2 \< A22 

i • (All + A22 + An + A22) < 1 
Vxi e SR 



Finally choosing Aii 



1 and A 



11 



1 \ * 
4 ' ^22 



i, this leads: 




4 ■ a;2 + 2 • a;i • X2 1< 7 



The plot for this region is depicted in Figure [TJ On the other hand, using 
Corollary [71 the domain of attraction is presented in Figure [51 

Example 2 (The Van der Pol System) . The Van der Pol system is a very well known 
system which serves as a benchmark in many areas of science, in particular it is 
stable for negative values of the parameter fi: 

\ Xi (t) = X2 

I X2it) = -Xi + ^■X2-il-xl) 
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Figure 2 . Domain of Attraction Comparison using Corollary [71 



Figure 3. Region for Hurwitz eigenvalues in the Van der Pol system. 



-''^iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiMiiiM 


iiSfss, 













Figure 4. Domain of Attraction obtained with Corollary [5j 



For this example, the systematic Corollary [7] will be considered. To apply this 
result, it is necessary to obtain the jacobian: 



df{x) 



1 

2^ ■ Xi ■ X2 /i • (1 - 



dx 

The eigenvalues are Hurwitz in the region depicted in Figure [31 It is clear that 
the stable region for the eigenvalues contains the equilibrium at the origin. In this 
way, the zone where the eigenvalue Xr in Corollary [7] is less than one is plotted in 
Figure m for /i = —1. 

Since this system is siii/jO the obtained stability region is smaller than the whole 
domain of attraction depicted in Figure [S] 

Example 3 (The Markus-Yamabe conjecture, discussion). One special class of ODE's 
are those with Hurwitz jacobian matrix for all 3?". For these systems, the intuition 



'^in the sense that possesses two main time-scales with different magnitude producing a sudden 
change into the direction motion of the system-this is one reason why it is used as a benchmark 
for numerical ODE solver testing 
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FiGURE 5. True Domain of Attraction for the Stable Van der Pol system. 

could say they are always stable. However, as it is well known since 1995 this is 
not true for dimensions bigger than 2 in general cases (see (31] and [7]). 

Taking a look of Corollary [71 is clear that the region for stability is the intersec- 
tion of two main regions: 

• The region where the eigenvalues of are Hurwitz. 

• The region where A^j < 1 . 

A conclusion using the theory in this paper is that the Markus-Yamabe conjec- 
ture was uncomplete, missing the second condition above. 

Example 4 (An example where the test is not applicable). This example is due to 
Krasovskii (see [3|, pp. 32, [24| and [2|). This system has an equilibrium at the 
origin but the jacobian is unstable for 5R^: 

Xi{t) = X2 

X2(t) = -Xi + X2 ■ {xl + xl)^ ■ S«"-( (^5^^)^ 

Since the jacobian is unstable, it is not possible to find an stable matrix F as 
required by Theorem otherwise, in the view of Lemma [2] the conclusion would 
be that the jacobian is stable leading a contradiction. 

Example 5 (Domain of attraction discussion). Even when the focus on this paper 
is to derive an analytic tool to determine stability /attraction for Nonlinear ODE's, 
this example intends to discuss how precise is our prediction of the truly domain 
of attraction, in fact, the Normal- Hopf example in Section [3] is very conservative 
but it is expected that the general scenario is of acceptable conservatism. In fact, 
it will be shown that for many systems the prediction exhibits high precision. 

A class of systems with known domain of attraction are planar systems with 
circular unstable limit cycles. Those dynamics can be written in polar coordinates 
as follows: 

r(t) = h{r) 

e(t) = 1 

here h{r) is negative for a finite interval on r containing the origin (r = 0), this 
idea is depicted in Figure [5] 

Notice that in order to have f{x = 0) = 0, it will be needed /i(0) = 0. Moreover 
for a finite circular domain of attraction, h{ii) = with fi a finite positive value. 
In this way and calculating the Jacobian: 
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Figure 6. Shape of h{r) 



dfjx) 
dx 



r ■ h{r) + cos{9Y ■ {h'{r) ■ r - 
r + siniO) ■ cos{6) ■ {h'(r) ■ r - 

_ dh(r) 



h{r)) 
-h{r)) 



-r + sin{9) ■ cos{9) ■ {h'(r) ■ r 
r ■ h{r) + sin{9Y ■ {h'{r) ■ r — 



-hir))- 
h{r)) 



where h'{r) = ' . Since the border of the domam of attraction is circular, in 
order to analyze the condition in Corollary [71 Aj^ < 1, it can be settle: 

r = ve* e K 

In this way, the Jacobian for 6* = (3: 



df{x) 
dx 



h'{ti) -I 
1 



then: 



Moreover: 



Finally: 



P 



1 1 
WJjTj 2 

1 2+h'ip.)'' 

2 2-/i'(/i) 



R = 2-n-{\P 



"0 


0" 




'1 


1" 


1 


1 


+ 









^1) 



A, 



Calculating Xb. < 1 it can be seen that for the interval /i G [0,0.2321], the 
estimation is perfect. Realize that this result is valid for any function h{r) making 
valid the result for a wide range of systems. 

5. Conclusions 

Sufficient conditions for asymptotic stability and attraction for nonlinear ODE's 
were presented. These conditions are systematic to apply, moreover the main the- 
orem provides a region where the ODE under analysis is asymptotically stable and 
a quadratic Lyapunov function. Realize that such a region serves as an estimation, 
sometimes coincident, of the domain of attraction and the Lyapunov function can 
extend this region by applying methods like the one in [6]. 



^The rest of the cases lead equivalent results since is a circular domain. 
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On the other hand, while Theorem |4] needs the input of a stable matrix F, 
Corollary [7] utilizes the jacobian matrix of the given nonlinear ODE. The former 
result provides a subset of the domain of attraction but the estimation depends 
strongly on the choice of that matrix F. In this sense, the corollary was presented 
to overcome this inconvenient making the method systematic. 

Several examples were presented for both Theorem |4] and Corollary [71 showing 
the applicability of such a results, moreover Example 5 analyzed systems where the 
estimation of the domain of attraction is exact. 

Many future directions can be depicted, in fact the eigenvalue Xr{x) in Corollary 
[7]it seems to increase from x — until it reaches the value 1 exactly when x belong to 
the border of the subset of the domain of attraction obtained. If this is confirmed, 
then this quantity A_r serves to determine the borders of the estimation of the 
domain of attraction when is equal to one. 

Besides, the theory supporting all the development in this paper is based on 
Theorem [1] which in essence is a sufficient condition for attraction in CPWL ODE's. 
Notice that the theorem is proved utilizing a finite improper integral, the finiteness 
is proved by requiring that all the matrices in each simplex of the CPWL under 
analysis being stable. This requirement could be in principle a bit conservative and 
should be improved allowing some of the matrices being unstable, then enlarging 
the attraction region provided by Xr. 

It will a topic for future research also the application of the methodologies on 
this paper to stabilize control systems also known as non-autonomous. In the case 
of control systems, one has to derive a controller but also to ensure that the system 
is stable at least for a finite region, in this sense, sufficient conditions for stability 
are the perfect scenario in order to obtain such a controllers. 

6. Appendix 

We present here the proof of those lemmas and theorems along the paper. 
proof of Lemma \^ Starting from the condition: 

(6.1) \-F + A + B\<X* 

In this way and considering a matrix P' = [Pi, P2, . . . , Pn]' , Pj — [Pji,Pj2, ■ • ■ , Pjn] 
and multiplying both sides of equation ()6.ip by | Pkj |, j, fc = 1, .., n: 

i\Pkj\-\ -F,i+A,i+B,i |<|Pfc, I -Tji 

In matrix form: 

(6.2) \ P \' ■ \ -F + A + B \<\ P \' -X* 
In a similar way is also true: 

(6.3) I i-F + A« + P«)' I • I P |< A*'- I P I 
Summing equations (16. 2p and (|6.3p : 

\P\'-\-F + A + B\ + \-F + A + B'\-\P \<X*'- I P I + I P r -A* 
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Realizing that | P ■ {-F + A + B) + {-F + A + B)' ■ | P |<| P 1' • | {-F + A + B) \ 
+ I (-F + A + By\-\P\ and muhiplying by | a; |, x £ at right and left: 

I X I' •(! Pi-F+A+B)+{-F+A+By-P 1). \x\<\x |' -(A^'- \ P \ + \ P \' -A^)- | x 
Taking into account that a;' | (•) | x <| a; |'| (•) || x |, it is possible to write: 

x'- I {P' ■ i-F + A + B) + {-F + A + By-P)\-x< 

X I' -(F'- I P I + I P r -A^)- I X 

Finally resorting to the property |a + 6|>|a| — Va, 6 £ 3? and summing 

x' ■ [F' ■ P + P' ■ F) ■ x: 

(6.4) 

x' -{P' ■A + A-P)-x < 
X I' -(F'- I P I + I P I' -A^)- I a; I +a;' • (P' • P + P' • P) ■ x) + x'- I (P' ■ P + P ■ P) I -a; 

where the property | a + 6 |> a + b,\fa,b £ 5R was invoked. Since P is stable 
if and only if there exists P = P' > such that {P' ■ F + F' ■ P) = -I, with / 
the identity matrix (see [25j, pp. 22), also considering the definition of positive 
definiteness (see [3]), in order to guarantee that the matrix A is negative definite 
(or stable) is enough to require: 

(6.5) _ _ 

I X I' -(A*'- I P I + I P r -A*)- I X I +x'-{F'-P+P'-F)-x)+x'- 1 P' B + B P \ -x < 
Recalhng that F' ■ P + P' ■ F = -I, then equation (|6?5|) looks like: 

{ |a;|'-(~'-|P| + |P|'-~)-|a:|+a;'-|P'-_B + B-P|-a; ^ ^ 
x' -x — 
x' • X > 

Utilizing the condition: 

(6.6) |P|<A* 
The requirement for stability of A will be: 



i.r. aA*+A*)'-iPi + ipr-(A*+A*)) .i.i 

(6.7) { - < 1 

' x'-x — 

x' • X > 

Since this inequality should be satisfied for all x £ 3?" in virtue of the definition 
of positiveness, the worst case leads: 



fa «^ \x\'-R-\x\ 

(6.8) max^ , „ , r ^ ^ 

I X I' • I X 

This is called Rayleigh's quotient and have important properties (see Jt\). To 
finish, realize that matrices R and | P | are positive (in the sense that all the 
elements are positive), so posses positive eigenvalues with positive corresponding 
eigenvector (see [3], PP- 72, Theorem 18). 
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Figure 7. Scheme of the Region Containing the CPWL approximation. 
In this way, the condition ensuring stabihty for A leads: 

Afl < 1 

This completes the proof. 

□ 

proof of Corollary Under the hypothesis of Lemma [2l and recalling the proof of 
this Lemma, in particular equation 16.41 leads: 

x' -{P-A + A' ■P)-x <0 

This completes the proof. 

□ 

proof of Theorem [7) The proof in what follows starts with the CPWL approxima- 
tion bounds introduced in j 23] : 

|/(x) - A« -x + bWI < A 

where A^*) -x + B^^^ are the linear approximations for /(x) in every simplex. This 
means that A^^^ ■ x + B^^^ is bounded by two vector fields: A + f{x) and —A + f{x), 
Figure [7] is depicting this idea. Notice that the dash region in that figure is the 
area containing the possible CPWL vector fields for a given A. 

Since the focus is to create conditions over f{x) related to {A*, A*}, we can 
consider the "worst case" A^~'^'> ■ x + B^~-^'> and A'-^^ ■ x + B'^^'^ and show that those 
planes define a convex set@ as follows (see Figure |S]): 

(6.9) A*"-^) • X + B^-^^ < f{x) - A < • X + B« < fix) + A < A^^^ ■ x + B^^ 

The key idea it is to show that the inequalities above, are satisfied for at least 
one point inside each simplex and centering the attention in a generic simplex with 
vertices: l^sJi"', z = l,...,n-|-l. In this way, if the gride size tends to zercP^. 
it is possible to prove that: A^'^ , B^"^ . 

Then, in place of using matrices B*^*^}, one uses the pair {A'^^^\b'^^^^}, 

which can be obtained for any generic nonlinear vector field f{x) without the ne- 
cessity of the evaluation of any CPWL approximation which moreover has to be 



^See [39|, pp. 160-165 
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\ 

• A™.x+B'"" 



Figure 8. "Worst case" planes. 



reduced more and more in order to consider the limiting case when the grid size 
tend to zero. 

Once the closed-form expressions for the pairs have been ob- 

tained, it is possible to analyze attraction or stability as it going to be shown in 
what follows. For the sake of clarity we will only consider the analysis for a; > 0, 
the case a; < shares conclusions. 

analysis: In this way, let's calculate first these worst case planes, then based 
on Figure m 



^i = /(^i) + A, P2-/(T4)-A, P3 = /(F2) + A, Pi^f{V2)-X 
f{Vi) + A = A(^) • Vi + S(^) 
f{V2) + A = A(^) • V2 + S(^) 
/(Vi) - A = A-(^) • Vi + S-(^) 

In this way, the matrices are given by: 



(6.10) 



(AW 



f{Vi)-f{V2) 
^-(A) ^ f{Vi)-f{V2) 

5-(A)^_[)ykp)].y,+/(^,)_A 

At this point, the following conditions need to be verified: 



(6.11) 



^(-^) • a; + S(-^) < fix) - A 
f[x) + \< A(^) -x + B^ 
> 



Using the notation Vi 
above yields: 



V2 = <5j; < without loss of generality, the conditions 
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• X + < fix) - A 4=» fix) - fiVi) < A^^'> ■ ix - Vi)) 

fix) + A < A(^) • X + S(^) ^ fix) - fiVi) > • ix - Vi) 

yx > 

Considering equation (|6.10p . it is clear that A^^^^ = A^~-'^\ then: 



< 2 • A 

fix)-f{Vi)<AW-ix-Vi)) 
fix) - fiVi) > ^(^) ■ ix - Vi) 

yx > 

First condition is trivially satisfied taking into account that a CPWL approxi- 
mation, always work with A > 0, however, second condition means: 



f fix) - fiV,) = AW ■ ix - Fi)) ^ ^gr^ff) = AW 
[Vx > 

Considering the limiting case Vi — J> V2, then: 

{ IMtJSyil - lim.r „ - Mill „ 

\ix > 

In the view of the mean value theorem, this condition is always satisfied if fix) 
is continuous and with first derivative continuous - see [20] . pp. 148. 

Once the validity of inequality ()6.9p has been proved for at least one point x £ 
Wi, V2], it is missing to prove the following: 

• If -)> V2, then: ^ A^*), ^ B^'). 

• The pair satisfies: | - F + A^^^^ + B^^^^ |< A*, | |< 
A*. 



The first item is going to be analyzed using inequality (|6.9p in the vertices 
{Vl,F2} and recalling that A^ = A(~^): 



using absolute values: 



< 



Vi 1 
V2 1 



< 



Vi 1 

1 



T ■ 

where T — 



AW - 
h - 



4^ 



<| T-^ I -T- 



^1 1 
V2 1 



and /i = -AW . Vi + /(Vi). In this way, if A ^ 0, then: 



AW - 

h - 



= 



in other words: A^ A^''> y B^"^ . 
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The second item above, it going to be proved in what follows and states: 

f| + |< a* 

jl |< A* 

For the sake of clarity in the exposition, it will be analyzed separately B^^^'' and 
For B'^^^^ : The requirement implies to satisfy the following: 



or, in other words: 



fiVi)>X+[^^^]-V,-X*, 



(6.12) 

For -F + A^"^^^ + B^*-^); This condition requires: 



• < _^ + ■ (1 - Vi) + f{V,) X 

-F+£^^^^-il-V,) + f{V,) + X<Xl 



This yields 
(6.13) 



ifiVi)<Xl+F-l^^:^^-{l~V,)~X 
\f{Vi)>~Xl + F-l^^^l^-{l-V,) + X 

Finally a region where the analysis is valid, has to be obtained. 



Validity Regions in 3?^: Equations (|6.12p and (|6.13l) provide a region where 
the analysis is valid or in other words where there exists a CPWL approximation 
satisfying: 

f| -F + |< A* 
jl |< A* 

In this way, the conclusions for stability and attraction valid for the CPWL vector 
field: A'-'^ ■ x + B^'i are also valid for A^^-^^ ■ x + S'^-^) which can be immediately 
translated into conclusions for the nonlinear ODE. 

For this reason, it is convenient to write together equations (16.121) and (|6.13p to 
give: 

/(y^)<_A + [mi5^].yi + A* 

f{v,)>x + [i^y^^^]-v,-xi 

f{Vi) <Xl + fix = 1)- ■{l-V,)-X 
f{V,) > -XI + fix ^ 1) - ) • (1 - Vi) + A 

The analysis will focus two main conditions: 

• The bounds for /(Vi) contain equilibrium points. 



(6.14) 
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• Suitable values for A in order to assure that a CPWL exists for our valid 
region in order to validate the proof using the CPWL theory theorem de- 
veloped in Section [T] 

The first item means a positive superior bound and a negative inferior one, per- 
forming this analysis to the conditions in (j6.14p . it will be also obtained appropriate 
values for A: 



Since A > 0, considering the limiting case Vi — > V2 and naming Vi = x, this 
leads: 



ST' - 
X . 2m < A* 

(l-x).^<A|+/(.x=l) 
.il-x).^>~Xt + f{x^l) 

In summary the region where the analysis is valid can be recast as: 
/l-^l<A2 

\|(l-x).^-/(:z: = l)|<A^ 
This completes the proof for the scalar case. 



5R" analysis: The development in this section extends the previous result for 3?^ 
to vector fields in K", to do this, we consider 2 • n "worst case" planes defining a 
convex set as follows: 



^A^r^) . X + B^r^^ < f,{x) - A, < • X + B^p < f,{x) + X, < A^^ ■ x + fif ^ 
J = l,2,...,n 

As for the 3?^ case we need to show: 



(A^f'^-x + Bp'^<A^P^-x + Bp^ 
A<^-^^ -x + Bp^^ <f,{x)-X 
(6 15) /,(-) + A<Af + 

i = l,2,...,n 

X = [xi,X2, . ■ .,X„]' 

yx, > 

In this way, it is first needed ^j-*'^' and B^^^ which can be obtained using n+l 
points {xi,X2, • ■ • , Xn+i} as follows: 
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" x'l 1" 




G K 










i = 1, 


. . . ,n 











(±A)'' 
/(±A) 



fc = 1, . . . ,n + 1 

The set of points {yi, y2, ■ ■ ■ , yn+i} is related to {xi,X2, • ■ • , x„+i} by: 

A = [Ai, A2, . . . , A„]' 
J ^l,...,n + l 

where {Vi, V2, . . . , V'n+i} is the set of vertices containing the general simplex i 
considered for the analysis. 

Then, it is possible to explicitly write A^^ and B^^ in matrix form: 



th 



(6.16) 



A 



(±A)' 



(±A) 



B 



" VI 1 


-1 







/.(^i)±A, 
/,(K+i)±A, 



Next, in order to effectively calculate A^^^'^ and Bj^^' from equation (|6.16p we 

■ VI 11 
need the inverse of the matrix : 

A useful way of doing this is the one presented in |17) , pp. 97, where applying 
row operations as in the well known method of Gauss for solving linear systems of 
algebraic equations it is possible to transform a matrix M in AI as follows: 



j{±A) 



M 




M 




where / is the identity matrix. Then by subtracting the first with the second 
row in M , the second with the third and so on, we arrive to: 



(6.17) 



M 



( VI - vi 0\ 
vi - Vi 



1-10 
01-10 
00 1-10 

00 ... 

00 



... 

... 

... 

1 -1 

1 



On the other hand is possible to choose n vertices forming a basis (see [23^ for a 
further reading): 
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14+1 = [■^i V2 



then there exists 7* G 5ft" such that: 



71 
72 

17" 



-K+i - [V^i - V^2 V2-V^ ... K - K+i] ■ 7* 
Now according the boundary confiauratiorP^ in [23 , pp. 58-65, equation (3.8) 
and more general in section (3.4.1), then it is possible to write every vertex for any 
simplex as follows: 

\i — I, . . . ,n 

where Cj G 0, 1, rrii = 1, 2, . . . , iV^ with Ni the number of simplicies used in each 
coordinate and Sx^ their grid size . In this way: 

'V.-V,+,= 4. ■ (E;=i(-1)^^ ■ - E ■=i(-l)^^ ■ - (-1) • e,+i) 



providing the present analysis for > 0, i = 1, 2, . . . , n-the rest of the cases are 
of similar consideration. This fact yields: 



'Sx^ 
<5.., 












In this way the matrix M from equation (|6.17p becomes 



M 



\r' -Q' 1 / 

Performing row operations we have 



1-10 
1-10 
1-10 








... 

... 

... 

1 -1 

1 



M 



( 



Vo 





1/ 



\ 

/* : 


-1 

V-7*'-/* 1+7:7 



^ ^This is how the set of vertices is called in 1231 
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where: 



1 -1 
1 -1 





Realizing that -7* • /* = -[71* (-71* + 72*) (-72* + 73) ■ • ■ (-Jn-i + ln)].^e 
obtain the desired inverse: 



V{ 1 



T- 

T- 



. . . 0' 

... 

••• it 

1 



Finally, A^^^'' and arise: 



1 -1 
1 



-71* -(-71* +72) 



... 

-1 ... 

1 -1 

... -(-7;.-i + 7,*j (1+7,:). 



(6.18) 



A 



i±\y 



/,(Vi)-/,(V2) 

/,(V^2)-/,(V3) 



* — with V(„_|_i)j each component of the vertex n + 1. Once 

and Bf^^ are already determined, the conditions for a convex set depicted 



where 7* 



in equation (|6.15l) lead: 



■a\-'^ ■ X + B^r^^ < ) . X + Bf^ ^ < [Bf^ - B-^'^) ^ 

2 • Aj > 
/,(x)-/,(K+i) = (a;-K+i)'-A* 
J = l,2,...,n 
As for 5R^, considering the limiting case Vi — ^ X^+i, then: 



A 



^-^^ ■ X + B^r^^ < ) • X + B)^' ^ < (B)^' - Bp^>) ^ 



f,{x) - /,(K+i) ^{x- K+i)' ■ ^U=v„,i 

The first condition above is always satisfied since all the CPWL approximations 
used along the proof provides a positive quantity Xj, on the other hand the last 
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requiremnet is also true for continuous vector fields f{x) and with continuous first 
derivative in the view of the mean value theorem for vector fields in 3?" (see [20] , 
pp. 148). 

On the other hand as for 3?^, the following it is remaining to be proved: 
• lfViV2, i= l,...,n, then: -^ylW, B^-^^l ^ B^'K 

l(.±A) p(±A) 



{A(±^),i3(±^)} satisfy: | F,k^A^f^^-B]^^> \< A*, 



l,...,n. 

The first item can be proved utilizing the inequality (j6.15p in the n ■ 
Vi, i = 1, . . . , n + 1 and noticing that A^^^ ^ A^"^) (see ((6?T8| ): 



1 vertices: 



A 



(A)' 



using absolute values: 





-Afy ~ A<f' 




T- 








. - Bf . 





where: 



< 



r v{ 1] 








\ ^1 ^1 


Vi 1 




-Afy 


< 


Vi 1 


.K+i 1. 








Yn+l 1. 



< T- 



h., - B^ 



(i)' 




In this way, if A — > 0, then: 



1 



A 



< |T~M ■ T ■ 



4'^ - Y 
h,~Bf 

in other words: v4(^) A^*) y S^'). 
The second item yields: 



ii)' 



A 



B 



(i) 
jk 



B^^ \< A 



•jk 



Jk 



J,k = l,...,n 

As for 5?^, the analysis will be split into and Fjk — A 



For b\^^^ . In this case, the following should be satisfied: 



(±A) 
jk 



B 



(±A). 



/j(K+i) + A,-K+i-^* <a;^ 

^/j(K+i)-A,-K+i.A* >-A* 



Realize that l^+i • A* play the role of 



f{Vi)-f(V2) 



for the case. 
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In this way: 
(6.19) 



For fj{ek) — ^jk^^ ^ ^j"^^^- guaranteed: 

k — 1, 2, . . . , n 
Recalling equation (|6.18p : 



A*k + lAVn+i) + A, - • A* - f,{e\) < A*, 
A*k + hiVn+i) - A, - K+i • A* - f,{e\) > -A*, 
k = 1,2,. 



this yields: 



'fjiVn+i) < \*k + - A, - A*, + K+1 • ^* 
(6.20) \ f, (K+i) > -A*, + /, (el) + A, - A*, + K+i • A* 

k = 1,2, ... ,n 

Validity Regions in K". Comparing equations (j6.19p with (j6.12p and (j6.20p with 
(|6.13l) for 5R" and respectively, we immediately see that the conclusions for 3?^ 
can be readily extended to the present 5R" analysis to give: 



V(n+l)i > ^ 

/,(14+i)<-A,+K+i-^; + A*fe 
.6 21) l/,(K+i)>A,+K+i-A*-A*, 

]fj iVn+i) < A*, + /, (efe) - A, - A*, + • A* 
fjiVn+i) > -A*, + /,(efc) + A, - A*, + • A* 
.j,^ = l,2,...,n 

As for K^, the present SR" needs to fit the requirement of the superior bound 
being positive and the inferior one negative in order to contain equilibriums of 
f{x). In this way: 

'A-fc-A,+K+i-A*>0 

+ A, + K+i ■A*<0 
< + f,{el) - A, - A*, + K+i ■A*>0 
-A*, + /,(efe) + A, - A*, + • A* < 
J,k^l,2,...,n 

The appropriate values for Aj arises: 
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A, < a;, - • A* 

< A, <A;, + /,(e,)-A*, + K+i-^* 
A, < A*,~/,(e,) + A*,-K+i-^; 
= 1,2, . . . ,n 

Since \j > 0, writing in matrix form, considering the limiting case Vi — )■ Vi+i 
and naming in general Vn+i — x, it is obtained: 



„*' dfj (x) ^ _ \* 
dx - _ -^j 

„*' dfjjx) < \* 
dx - ■^J 

dfjjx) 



[-X* + I) ■ ^ > -A* + Fk 
i-x*' + I) ■ ^ < A* + Fk 
j ^l,2,...,n 



where ']' and F — [Fi F2 ...F„]. The conditions above 

can be compactly written as follows choosing the minimum bound for the first two 
inequalities above: 



(6.22) 



9 fix}' 



dx 



< minj=i^2,...,n A* 



df(x)' _ 



r dfijx) 

dxi 
dh(x) 



dx2 



dfijx) 
dx 



F\<X 

dh(x) 

dxi 
df2(x) 

dX2 



dh(x) 
dx 



dfrAxY] 

dxi 
8x2 



dU(x) 



These conditions can be further simplified as follows: 



(a OQ^ / I • I ^ '^«'^J = l,2,...,n A* 



To complete the proof, it is remarkable that the conditions obtained in this proof 
are independent of the coordinate system chosen. 

Up to this point, it has been proved that there exists a region given by (j6.23p and 
a precision CPWL under which all the matrices are Hurwitz in each simplex, 
because of that and since the objective it is to prove that the aforementioned region 
it is also a region where the nonlinear ODE: x{t)f{x) is attractive to its equilibrium 
point, it will be necessary to prove that this region is an invariant set. 

By virtue of inequality (|6.15p . if the trajectories of: x{t) = A^*) • x + _B^*', are 
considered: 
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" v ' 

x(t) 

I B(±^) |< mm,=i,2,...,„ A* 

= ^•x|<min,=i,2,...,„ A*, | ^ _ F' |< (A* ~ A*)} 



f(Vi)-f(V2) 

/(V^2)-/(V3) 

5xo 



<5x„ 



In order to solve the differential inequality above, the technique in [15] is going 
to be applied to give: 

' min{i_2} {x^*{t),x'^*{t)} < x{t) < max{i_2} {x^* {t) , x"^* (t)} 
x^*{t) = A(^) • x^*{t) + mm,=i,2,...,„ A* 



„2* 



(i) = A(^) • x^*{t) - mm,=i,2,...^„ A* 



It is possible then, to assure that the trajectories given by the inequality: A^-^^^ ■ 

X + < A(^) • X + B^"^ < A(^) • X + B(^) are contained in VL, "it > Q if it is 
. ' ~ 

x(t) 

ensured that the bounds are also contained there. 

To do this, the extremums of x^''^*{t) are looked for, that means, the instants 
of time for which x^''^*(t) reach their maximums and minimums. In order to find 
these maximums or minimums, the elemental calculus result for which the instant of 
time where the extremum occur are utilized, in this way, it is necessary to evaluate 
x^''^*{t) = and does not exist. 

It is clear that the first possibility only happens if t — 00 in the view that A(^) 
is Hurwitz, while the second one indicates t = due to the discontinuity introduced 
in considering only positive evolutions in time: {t G 

These reasons and considering the initial condition x^'^*(0) always inside the set 



il, it will be enough to prove that: lirrit- 



*{t) £ fl, that is: 




equivalently: 

In this way, if A 0: 



x^*{t) = -mmi^i,: 

x'^*{t) = • 771111,^1,2,. 



yl(^) • x^*{t) = -A* < mm,=i,2,...,„ A* 
^(^)a;2*(i) = A* < TOm,=i,2,...,„ A* 



df{x) 



dx 



x \< min,=i,2,....,n A* 



This completes the proof. 



□ 
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proof of Lemma[^ As it was proved in Theorem |4] each matrix A^^^ is Hurwitz. 
In this way and applying the Corollary 4.1 in [22] it is immediately seen that 
V{x) = x' ■ P ■ X is a common quadratic Lyapunov function of: 

xW(t) ■x^''> +B^'\ \/i^l,2,...N 

where N is the number of total simplices used for the partition of the CPWL 
ode's considered. If, on the other hand, §7 = and using the bounds introduced 
in [23], it is possible to write: 

I /(xW) - . _ |< A, Vi = 1, 2, . . . iV 
Multiplying both right sides by | a;*^*^ | • | P |: 

I I • I P I • I /(a;(*)) - A« • a;(*) - |<| a;^' | • | P | A, Vi = 1, 2, . . . TV 

Since | x«' • (P • /(a;«) - A« • a;W - sW) |<| a;W' | • 1 P 1 • | /(xW) - • 
_ 1^ it is obtained: 

(6.24) I x^^y - P- (/(xW) - -2;^ -P(*)) |<| xW I • I P I A, = 1,2,...A^ 
With a similar procedure is also valid: 

(6.25) I (/(a;(*)) -AW ■ - P«) • P ■ x^^^ |<| x^'^' | ■ | P | A, Vi = l,2,...iV 
Finally summing equations (|6.24p and ()6.25p : 

2-a;(*)'-P-/(x(*))-2-x(^)'-P-(AW-x(*)+pW) < 2- | x^'^ | • | P | A, Vi 1, 2, . . . 

Since V = x^'^ ■ P • x^-^^ is a Lyapunov function for ^^'^ • x^'^ + P*^*', then 
2 • • P • (A(^) • + P(')) < 0. This is showing that: 

2-a;W-P-/(a;W)-2- | x^*)' | • | P | A < 2-2;(*)'-P-(^(*)-x(^)+p(^)) < 0, Vi = l,2,...iV 

If the limiting case A is considered, then in virtue of Error Bounds the 
Dynamics introduced in |16] . it is known that the trajectories x'*^ x{t), where 
x{t) — f{x), this leads: 

limx^o 2-x(^)'-P-/(x(*))-2- | x'^'^ | • | P | A < ^ 2-x{t)-P-f{x{t)) < 0, Vi = 1, 2, . 
This completes the proof. 

□ 

proof of Corollary As for the proof of Theorem [H it will be needed a Hurwitz 
matrix F such that: | - P |< A* - A*. In this way if is Hurwitz for all 

X e il and using A* = A*, then if ft is as follows: 
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'n = {x:XR<l} 

R=2■{\p\■\^^■x\■[l 1 
- 2^ ' 

H(x) 

with the biggest eigenvalue of R 

the ODE: i{t) = fix) is attractive to the origin (equihbrium point) in the view 
of Theorem |4l 

This completes the proof. 

□ 
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